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ABSTRACT 

A detailed analysis of the optical polarimetric variability of the TeV blazar lES 
1959+650 from 2007 October 18 to 2011 May 5 is presented. The source showed 
a maximum and minimum brightness states in the R-band of 14.08+0.03 mag and 
15.20+0.03 mag, respectively, with a maximum variation of 1.12 mag, and also a 
maximum polarization degree of P =(12.2+0.7)%, with a maximum variation of 
10.7%. From August to November 2009, a correlation between the optical 7?-band 
flux and the degree of linear polarization was found, with a correlation coefficient 
'"po/=0.984±0.025. The source presented a preferential position angle of optical polar- 
ization of ~153°, with variations of 10°-50°, that is in agreement with the projected 
position angle of the parsec scale jet found at 43 GHz. From the Stokes parameters we 
infer the existence of two optically-thin synchrotron components that contribute to the 
polarized flux. One of them is stable, with a constant polarization degree of 4%. As- 
suming a stationary shock for the variable component, we estimated some parameters 
associated with the physics of the relativistic jet: the magnetic field, B ~0.06 G, the 
Doppler factor, 6q ~23, the viewing angle, O ~2.4°, and the size of the emission region 
riy ~5.6xlO'^ cm. Our study is consistent with the spine-sheath model to explain the 
polarimetric variability displayed by this source during our monitoring. 

Subject headings: (galaxies:) BL Lacertae objects: individual (lES 1959+650) — galaxies: 
jets — galaxies: photometry — polarization 
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1. Introduction 



Blazars constitute the most extreme subclass of the active galactic nuclei (AGN). They 
are radio-loud AGN and include the BL Lacertae objects and the Flat Spectrum Radio Qua s ars 



(lAngel & Stockman 



Ghisellini et al 



19801 : 



Impey & Tapia 



19901 : 



Maraschi et al. 



19951: 



Fossati et al. 



19971 : 



19981). BL Lac objects show yery weak emission li nes with equivalent width of 



less than 5 A or no emission lines at all (e.g. 



Ghisellini et al 



19981). In contrast, flat- spectrum 



radio quasars (FSRQ) exhibit broad optical emis sion lines . This suggests that in BL Lac objects 



accretion disk emission is strongly suppressed (Levinson 



20061) . Blazars show strong flux 



variability, s uperluminal motion , and a non-thermal continuum extending f r om radio to TeV y-r ay 



regions (e.g 



Marscher et al 



Valtaoia et al 



201( 



1990; 



Abdo et al. 



Sikora & Made] ski 



2011 



2010bh. 



20011: 



Bottchei 



2007£ 



Agudo et al.l 



BOUallbl: 



In the spectral energy distribution (SED) of blazars, the non-thermal continuum emission 
shows two broad low and high energy-peaked components. It is widely accepted that the 
low-energy peak component is produced by syn chrotron radiation from a relativistic jet ([Rees 



1967 



Jones et al. 



19741 : 



Marscher & Gear 



1985h . However, the nature of the high-energy peak 



component is more controversial. There are two different approaches to explain its origin. The 
first one comes from the so-called Leptonic models, which are based on the Inverse Co mpton 



(IC) scattering of soft photons by the same electrons emittin g the synchrotron radiation (|Konigl 



1981 



Marscher & Gear 



1985 



Ghisellini & Tavecchio 



20091). In these models, if the seed photons 



are provided by the synchrotron radiation emitted at lower energies b y the same IC scatteri ng 



electrons, the mode 



s are known as synchrotron self -Compton (SSC, 



Celotti & Ghisellini 



mm; 



Maraschi et al. 



( 1992 ) 



Marscher & Geaii (|l985|) '). If the dominant contribution to the seed 



photons field for IC comes from regions external to the jet, they are known as External Compton 



(EC) models. Possible sources of externa 



entering into the emission region directly (|Dermer et al 



seed photons include: (i) accretion disk p 



1992 



Dermer & Schlickeiser 



hotons 



1993|); (ii) 



photons reprocessed by the clouds within the broad line region (BLR) (|Sikora et al.lll994() : (iii) 
jet synchrotron emission r eprocessed by circumnucle ar material located close to the accretion 



disk and its X-ray corona (IGhisellini & Madau 



surrounding the central engine (JBlazejowski et al. 



1996h: (i v) infrared emission from a dusty torus 



20001) ■ The second approach used to explain 



the high-energy peak emission comes from the Ha dronic models wh ere it is assumed that the 



emission is produced by pairs and pion pro duction (iMannheim 



radiation from protons, n- and //- particles (|Aharonian 



2000; 



1993), as well as 



jy syn chrotron 



MUcke & Protheroe 



2001). In some 



cases, a hybrid model may be considere d when e 



the hadronic synchrotron mirror model (jBottcher 



ements of both models might be relevant such as 



2007bl) 



A distinctive feature of blazar emission is its high and variabl e linear polarization in radio 

d). 



and optical bands (lAngel & Stockman 



19801 : 



Impey & Tapia 



1990I) . This property associates 



the observed emission with beamed synchrotron radiation that is pr oduced by a relativistic 



jet viewed with a sniall angle relative to the observer's line of si ght (Bla ndford & Rees 



Blandford & Kdnigl 



1995 



1979 



Sikora & Madejski 



Blandford & Pavn 



2001 



Levinson 



982 



Urrv & PadovanJl995l : 



1978 



Padovani & Giommi 



20061). In recent years, th ere has been an increasec 



interest regarding the op t ical polarimetric properties of blazars (e.g. 



D'Ammando et al. 



Andruchow et al.l 



2009; 



2011 



D'Arcangelo et al. 



Ikejiri et al. 



2011 



20091 : 



Marscher et al. 



Hagen-Thom et al. 



2010 



Chandra et al. 



2008 



2011 



). This is due to the fact that polarization studies can 
provide useful information on the relative structure and configuration of the magnetic field 
associated with the relativistic jet. Variations of the position angle of the polarization vector 
may be associated to variations of the direction of the magnetic field's vector along the line of 
sight. Also, the degree of optical polarization could be related to the level of ordering of an 
initially tangled magnetic field or to the electron energy distribution within the emission region 



(|Angel & StockmanI 



1980 



Lister & Smith 



20001: 



Dulwich et al. 



20091) . 



In general, the linear polarization in blazars varies randomly although for some objects a 



systematic behavior of this quantity is observed during some period of time. The rotation of 
the polarization vect or might be produced by a magnetic field with a helical structure or by 



a bending of the jet (iSasada et al. 



2011 



Marscher et al. 



2010|) . It has been suggested that the 



polarization randorn behavior might be produced by several components of the polarized emission 



(Holmes et al. 



19841) . A succ essful rnodel us ed to explain the observed polarization in these 



objects has been proposed by 



Brindle 



(|l996l) . with the inclusion of two components. In such 



a model, the synchrotron sources of polarized radiation are optically thin and have different 



polarimetric characteristics. One of the polarized componen ts is stable wh i 



Kikuchi et al 



e the other shows 



(Il976l) . the authors 



a chaotic behavior. In a study done on the blazar OJ 287 by 

found evidence for two basic components: a steady polarized component, which they suggested is 

always present, and a randomly variable component, dominating during active phases. 

The TeV-blazar lES 1959+650 was discovered in the radi o band as part of a 4.85 GHz 



survey performed with the 91m NRAO Green Bank telescope (|Becker et al 



1991 



observations have been performed with the Very Long Baseline Array (VLB A, 



). Systematic 



Bondi et al. 



20011: 



Rector et al.ii2003i : iPiner & Edwards 



2005-2009 revealed spine-sheath str uctures in t 



fractional polarization distributions (|Piner et al. 



2004r) . Polarimetric and 43 GHz VLBA observations from 
le ele ctric vector position angle (EVPA) and 



2010). These authors concluded that the blazar 



lES 1959-1-650 consists of a compact core with a flux density of ~60 mJy and a ~1 milli-arcsecond 
jet extending to the southeast at a position angle of about 150°. On the other hand, the source was 



observed also in the optical bands where it displayed large and fast flux variations (ISchachter et al. 



1993 



Villata et al. 



2000; 



Krawczvnski et al. 



201o[ 



2004: 



Gutierrez et al. 



2006: 



Havashida et al 



2008 



Tagliaferri et al. 



2008 



Bottacini et al. 



Ikejiri et al 



2011 



Kapanadze & Janiashvilil 



20121) 



The source is known to be hosted by an elliptical galaxy at z = 0.04 7 and with Mr =-23. This host 



galaxy shows a disk and an absorption dust lane (IHeidt et al. 



19991). The SED of lES 1959+650 



shows its first synchrotron peak at UV-X-ray frequencies , therefore this object is classified as 



high-peaked BL Lac object (HBL) (IBottacini et al. 



20101) . The mass of the central black hole 
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has been estimated to be ~ 1.5 x IO^Mq (Fa 



been observed with RQSAT and BeppoSAXdBeckmann et al. 



XMM-Newton dGiebels et al 



2OO2I : 



omo et alJl2002h. In X -ravs. lES 1959+650 has 



Krawczvnski et al 



2004 



2002), with RXTE-ARGQS and 



Gutierrez et al. 



2006 ). These data 



showed that the synchrotron peak was in the range of 0.1-0.7 keV, and the overall optical and 
X-ray spectrum of up to 45 keV is due to syn chrotron emission with the peak moving towards 



high energy frequencies when the flux is high (JTagliaferri et al 



20031) . Also, the first y-ray signal 



at ver y high energies fr om lES 1959-1-650 was reported in 1998 by the Seven Telescope Array in 



Uta h (INishivamal 



(see 



Hartman et al 



1999h. Subsequently, the source was observed aga in emitting high energy y-rays 



1999 : 



Aharonian et al, 



2003 : 



Holder et al, 



20031) . In 2002 this object showed 



two T eV flares without simultaneous X-ray flares, a behavior sometimes referred to as orphan 



flares (Krawczvnski et al. 



2004 



Daniel et al. 



expected within the frame of SSC models (|Ghisellini et al. 



20051). Such orphan flares in VHE y-rays were not 



19981). After the detection of TeV 



emission, the source becarne a target for different multiwavelength campa i gns (.Krawczvnski et al 



20041 : 



Gutierrez et al. 



20061 



Tagliaferri et al. 



2008 



Havashida et al. 



2008 



Bottacini et al. 



2mm- 



This blazar is liste d in The First Catalog of AGN detected with Fermi Large Area Telescope 



(|Abdo et al 



2010a ) 



In this paper we report the results obtained from the photopolarimetric monitoring of 
the TeV-blazar lES 1959-1-650. Our main goal is to establish the long-term optical variability 
properties of the polarized emission in the R-band. Variability of the Stoke's parameters obtained 
from our observations is analyzed in terms of a two-component model. Estimations of some 
physical parameters that are known to be associated with the kinematics of the relativistic jet 
were obtained. Also, we present a comparison of the polarization properties found in our study 



Finer et al. 



( 2010 ). 



with recent radio maps obtained with the VLBI at 43-GHz on this source by 
Our work suggests that the observed radiation originates in a region inside the jet where a 
standing shock is produced. The paper is organized as follows: in section |2] a description of our 
observations and data reduction process is presented. In section [3] we show our observational 



results. In section |4] we present the polarimetric data analysis. In section [5] a discussion of our 
photometric and polarimetric variability results is presented and in section [6] we give a summary 
of our main results. Throughout this paper we use a standard cosmology with Hq = 71km s"^ 
Mpc-^ , Q,„ = 0.27, and Qa=0.73. 

2. Observations and Data Reduction 

The observations were carried out with the 0.84 m- f/15 Ritchey-Chretien telescope at the 
Observatorio Astronomico Nacional of San Pedro Martir (OAN-SPM) in Baja California, M exico. 



Sorcia et al. 



2011 



and 



We used the instrument Polima which is a direct image polarimeter (see 
references therein). All observations were taken as a part of a support program dedicated to 
the optical R-band monitoring of 37 blazars in the framework of the GLAST- AGI LE Support 

2 Polima 



Program (GASP) of the Whole Earth Blazar Telescope (WEBT) (e.g. iVillata 



20101) 



consists of a rotating Glan-Taylor prism driven by a stepper motor with an accuracy of 0. 1 ° . The 
Polima' s prism specifications show a transmittance of 90%, extinction ratio between the two 
states of polarization of 5x10"^, and a wavelength range of 215-2300 nm. Polima has a clean 
(unvignetted) field of view of (3x3)'in the sky plane. The telescope has an equatorial mount 
and its pointing angle spans from declinations between +70 and -45 degrees, and hour angles 
between -i-5 and -5 hours. We used three different CCD cameras during the monitored period 
with formats of (Ixl)k in two of them, and one with (2x2)k. Therefore, pixel sizes go from 
13.5 to 24/im and plate scales from 0.22 to 0.39 arc-sec pixel"^. Due to the small pixel size of 
the CCDs, and to decrease the integration and readout times, a binning mode of 2 x 2 was used 
in all the observing runs. The exposure time was 240 s per image for lES 1959-1-650. Polima 



'A detailed description of our photopolarimetric monitoring program on TeV Blazars can be 



found in http://www.astrossp.unam.mx/blazarel 
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is a single-beam device with a very slow modulation and quite sensitive to the sky noise level. 
The errors obtained with this instrument could be larger when the sky noise is high since it will 
dominate over the instrumental errors, so photometric conditions are required for an accurate 
polarimetry. The optical polarimetric monitoring of the blazar lES 1959+650 was carried out 
from 2007 October 18 JD(2454392) to 2011 May 5 JD(2455687). During this period of time, we 
carried out 25 observing runs of seven nights per run, centered on the new moon phase when the 
object was visible. In total we collected 106 data points. 

Four images with the 7?-band filter and relative position angles were taken. The sequence of 
0°, 90°, 45°, and 135° for the relative prism position angle was used to reduce the sky variations 
influence to calculate both linear polarization and flux from the object. Observations of nearby 
blank sky regions were obtained at different prism orientations using the same on-chip exposure 
times and on-source observations. Dark-current frames were not taken since the detectors were 
always operating at cryogenic temperatures. Flat fields were taken at the four polarizer positions 
in the R filter at dusk and dawn. Bias frames were also taken at the middle of the night. For each 
prism position, bias frames were subtracted from the flat field images. Then, flat field images were 
combined to obtain an average flat field image for each prism position. The bias frame was also 
subtracted from all object images and the resulting frames were multiplied by the mean value of 
the combined bias-corrected flats. Finally, the object image was divided by the combined flat field. 
Data were reduced with a pipeline especially written for our monitoring program (Hiriart 2013, 
in prep). Polarimetric calibrations were done using the polarized stand ard stars ViCygl2 a nd 



1992). The 



HD 155 197 and the unpolarized standard stars GD319 and BD-h332642 (iSchmidt et all 
instrumental polarization found is 0.6+0.5% and the zero position angle at (-90 + 2)°. Photometric 
i?-band magnitudes were determined from two orthogonal measurements: /i = /(0°) + /(90°) and 
/2 = /(45°) + /(1 35°) where f{x) is the flux of the object (or standard star) obtained at polarizer 
position X. The instrumental flux is the average of these two fluxes, from which instrumental 
magnitudes for each object were obtained. The magnitudes were measured using the aperture 



photometry technique. Then, we calculated the object's 7?-band magnitudes using comparison 
star #2, located ~2.1 'away from the targe t. The calibrated ma gnitude of the comparison star #2 in 



the 7? -band used is 12.53+0.02 mag, from 



Villata et al 



( 1998 ). 



We calculated the normalized Stokes parameters q and u for each object as: 



and 



u = 



/(0°)-/(90°) 
/(0°) + /(90°)' 



/(45°)-/(135°) 
/(45°) + /(135°) 



(1) 



(2) 



The q and u polarization parameters are related to the fractional polarization p and the 
position angle of polarization 6 by: 

p = ^q^ + u^ , 



and 



1 u 
6 = - arctan(-) 

2 q 



(3) 



(4) 



Defining the fractional po larization in this way, we are assuming that the circular polarization 



is negligible (see 



Jonesl(ll988h ) 



2.1. Error calculations 



From equations (H]) and (|2]), we propagate the errors for u and q to obtain: 






"(/45+/l35)2 



"(/45+/l35)2 



and 



/r ^-^90 -,2 , r 2/0 



(/0 + /90) 



'ifo+f9o) 



(5) 



(6) 
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where cTf is the error in the instrumental flux (see lNewberrylll99l|) and fx is the instrumental flux 
at the orientation x of the polarization analyzer. 



Similarly, the errors in the polarization degree, ep, and polarization position angle, ee, follow 
from equations (|3]) and dH), we obtain: 



£p 



v^ 



+ o-i 



(7) 



and 



^e 



= Y' V*^'?^")^ + (uo-gy , 



2p 



(8) 



where cTu and cr^ are given by eqs. (|5]) and dD, respectively, and p is the measured polarization. 



3. Observational Results 



3.1. Global variability properties 



The i?-band magnitudes were converted into apparent fluxes using the expression: 
Fobs = KoX 10-°^'"\ with Ko = 3.08 x 10^ mJy, for an effective wavelength of A = 640 nm. The 
photometry was done with an aperture of 3 " in all runs, and the obtained images have an average 
FWHM (Full Width and Half Maximum) ~3 ". The flu xes were corrected fo r the contribution of 
the host galaxy, subtracting 0.84+0.02 mJy according to 



Nilsson et al 



( 2007 ). 



We have to consider the ambiguity of 180° in the polarization angle. For this purpose we 
corrected the polarization angle assuming that the differences between the polarization angle of 
temporal adjacent data should be less than 90°. We defined this difference as: 



\m\ = K^l - e^l - V^(en+l)2 + (TiOnr, 



(9) 



where 6n+i and 0„ are the n + \ and n-th polarization angles and cr(9„+i) and cr{dn) their errors. If 



11 



\^e„\ < 90°, n o correction is need ed. If A0„ < -90°, we add 180° to 0„+i. If A^,, > 90°, we add 



-180° to 0„+i (ISasada et al 



201 ih . 



Therefore, all data corrected by the contribution of the host galaxy and the 180° ambiguity 
are presented in Tabled] There, we show in column (1) the Julian Date, in columns (2) and (3) 
the polarization degree and their errors, in columns (4) and (5) the polarization position angle and 
their errors, in columns (6) and (7) the i?-band magnitudes and their errors. The i?-band fluxes 
and their errors are given in columns (8) and (9). Figure [1] shows the 7? -band light curve, the 
polarization degree p and the position angle 6 obtained for all 25 observing runs in ~3.6yr. For 
clarity, the entire period of observations has been divided into three main cycles: Cycle I: from 
2008 May 04 to 2008 Dec 03 (JD 2454591-2454804). Cycle II: from 2009 Apr 22 to 2009 Nov 
18 (JD 2454944-2455154). Cycle III: from 2010 May 08 to 2011 May 05 (JD 2455325-2455687). 
These cycles are marked with dashed vertical lines in Figure [Hand will be discussed in more detail 
in the next paragraphs. 

To perform the statistical analysis on our data, nightly average values and the corresponding 
standard deviations were estimated. In order to estimate the variability in flux, polarization degree 
and polarization positio n angle, a y^-test was c arried out. The amplitude of the variations Y% 



was obtained following 

instead of magnitude diff"erences. Therefore 

100 



Heidt & Wagned (|l996h . but their relations were applied to flux densities 



y(%) = 



{S) 



^r 



'-' mm/ ^'^ c 



(10) 



where S max and 5„„„ are the maximum and minimum values of the flux density, respectively. {S) 
is the mean value, and al = crl^ax + '^mm- "^^^ variability is described by the fluctuation index // 
defined by 

fTc 

(11) 



"="»§* 



and the fractional variability index of the source T obtained from the individual nights: 



r 



S - S ■ 

'-' max '-' mm 

S* + S ■ 

'-' max ^ '-' mm 



(12) 
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Table [2] shows the results obtained from the statistical analysis. In column (1) the 
corresponding cycle is shown, in column (2) the variable parameters, from columns (3) to (10) 
we present the average value for each variable parameter, the maximum and minimum value 
observed, the maximum variation obtained b^max, the variability amplitude Y(%), the variability 
index //(%), the variability fraction T and the statistic ^^, respectively. 

We found for lES 1959+650 an average flux of (F) = 4.38+1.02 mJy, with a variation of 
Yp = 105.25%, corresponding to AF = 4.62 mJy; the fluctuation index is ixp = 23.28% and a 
fractional variability index is Tp = 0.47. The maximum and minimum brightness are 14.08+0.03 
mag (7.19+0.17 mJy) and 15.20+0.03 mag (2.57+0.06 mJy), respectively (considering all cycles). 

We have estimat ed the minimum flux variability timescale using the definition proposed by 



Burbidge et al. 



( 1974 : 



dt/\n(Fi/F2) 



(13) 



where dt is the time interval between flux measurements F^ and F2, with Fi > Fi- We have 
calculated all possible timescales Tij for any pair of observations for which | F, - F j \> o-f. + crp. 
at frequency v. The minimum timescale is obtained when: 



Tv = mm r, 



IJ,VI , 



(14) 



where i = \, ...,N - 1; j = z + I, ...,N, and A'^ is the number of observations. The uncertainties 
associated to Ty were obtained through the errors in the flux measurements. The minimum flux 



variability timesca 
in agreement with 



e obtained from ou r data in th e i?-band is ts =?,,„■„ =9. 7 4+ 1.1 7 d. This result is 



Krawczvnski et al. 



(2004) and 



Tagliaferri et al. 



(20081). 



Using the results of the observa tional runs presented here, we can quantify the duty cycle 



(DC) following 



Romero et al. 



( 1999 ): 



2"_, A^,-(1/Az',) 



(15) 
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where n is the number of nights, A^, is equal to or 1, depending on whether or not variations 
were detected in Ar, = Ati^obA^ + zy\ where Af, is the time interval at the source rest frame, and 
Ati^obs is the time interval measured between each run (at the observer's rest frame). We obtained a 
duty-cycle DC = 38.7% in the 7?-band flux for variations that are above 2cr, with cr = 0.27 mJy. 

3.2. Photometric variability 

In addition to our main three Cycles, we include in Table |2]the results obtained considering 
the entire data set. From these data, a maximum of R =14.08 mag is observed, and a variation of 
Amu =1.12 mag in At =329 d (or 0.9 yr) is found. Therefore, during our monitoring period the 
source shows a maximum scale brightness variation in time scales of A;^ ~ months. 

In Cycle I, the object showed a discrete variability with a relative flux maximum of 4.46+0.07 
mJy (14.60+0.02 mag) in 2008 May 04 (JD 2454591), and a relative minimum of 3.33+0.06 mJy 
(14.91+0.02 mag) in 2008 Dec 03 (JD 2454804). There is a moderate rise in flux between 2008 
June 06 (JD 2454624) and 2008 October 26 (JD 2454766) from 3.74+0.06 mJy to 4.32+0.06 mJy 
in 142 d, at a rate of 0.004 mJy/d. 

In Cycle II the object has been active for 148 d from 2009 April to 2009 September (JD 
2454944-2455092). The peak flux of the object was 7.19+0.17 mJy (14.08±0.02 mag) at 2009 
June 19 (JD 2455002). In this cycle we found a maximum brightness which is also a maximum 
value observed for the entire period of observations. At the onset of the activity, 2009 April 22 
(JD 2454944), the flux was 5.66+0.08 mJy (14.34+0.01 mag), then increased to 7.19+0.17 mJy in 
58 d (with a rising rate of 0.026 mJy/d). In 2009 September 17 (JD 2455092) it feU to 4.92+0.07 
mJy (14.49+0.01 mag) in about 90 d (with a decay rate of ~ 0.025 mJy/d). 

Finally, in Cycle III the source showed a constant increasing flux from 2.57+0.06 mJy 
(15.20+0.02 mag) in 2010 May 14 (JD 2455331) to 5.96+0.08 mJy (14.28±0.01 mag) in 2010 



-14- 
November 04 (JD 2455505) in 174 d (rising rate -0.020 mJy/d). 

3.3. Polarimetric variability 

3.3.1. Polarization degree variability 

To establish a possible correlation between the polarization degree and the i?-band fluxes, 
a Pearson's correlation coefficient was calculated {rpoi). This coefficient was tested through the 
Student's ^test. Using all data, we found that there is no correlation between the 7? -band flux 
and the polarization degree. The polarization degree showed a random variability behavior, and a 
maximum and a minimum of 12.2% (JD 2454656) and 1.5% (JD 2455685), respectively. So, the 
maximum variability observed was A/'=10.7 %, in A? =1029d (or 2.8 yr). 

From the light curve obtained for each cycle (Figure [D, the Cycle II was split into two 
periods that we called Cycle Ila from 2009 Apr 22 to 2009 Jun 25, and Cycle lib from 2009 Aug 
14 to 2009 Nov 18 (see Figure |2l). We cannot see any correlation between the i?-band brightness 
and polarization degree in Cycle Ila. However, for Cycle lib we found a strong correlation 
between these two parameters, with a Pearson's correlation coefficient of r^o/ = 0.984 + 0.025, see 
Figure |3] 

The maximum and minimum polarization degrees for each cycle are shown in the Table |2] In 
Cycle I, the maximum variability observed of the polarization degree is AP = 8.9% in A ? = 34d; 
in Cycle Ila, ^P = 3.5% in ^t = 30d; in Cycle lib, ^P = 5.0% in At = 90d; and in Cycle III, 
AP = 6.9% in A?=29d. 
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3.3.2. Position angle variability 

In general our data show that there is not a clear correlation between the polarization angle 
and the i?-band flux. However, the polarization angle shows a general trend to remain around an 
average value. We have looked for a possible correlation between the polarization degree and the 
position angle using the entire data set, see Figure |4] This plot shows that 9 has an average value 
of 153 + 16°, with maximum variations of ~ 50° around the average value. This suggests, that the 
position angle has a preferential position of ~ 150°, independently of the polarization degree. The 
preferred position angle of optical po larization is in go od agreement with the projected position 



angle of the parsec scale jet found by lPiner et al. 



(|2010h from 43 GHz images. We did not observe 



any significant rotation of the position angle during the entire period of observations. 

In Cycle I, the polarization angle rises from 135° (JD 2454594) to 181° (JD 2454681) at a 
rate of 0.53° per day, and falls to 132° (JD 2454707) with a rate of 1.88° per day; i.e. variations 
of 10°- 30° per run. When the magnitude falls to 14.75 mag (JD 2454773) the polarization angle 
rises to 186°. So, the polarization angle shows two oscillations when the magnitude increases. In 
Cycle II, the polarization angle is close to a constant value, in the period of maximum activity 
of the source. However, there is a slight change in the polarization angle from 146° to 164°, i.e. 
Ae -18°, when the magnitude faUs from 14.44 mag to 14.72 mag (JD 2455059-2455154). In 
Cycle III, we see a steady increase in brightness while the polarization angle 6 changes from 157° 
to 203°. It is worth to note that in this cycle, the amplitude variations of the polarization angle 
increase from 10° to 50° over the average value, as it does the brightness of the source. 



4. Polarimetric Analysis 

The study of the variability properties of the polarized emission in blazars allowed us to 
derive some properties associated with the variable source. In particular, for lES 1959-1-650 we 
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found from our observations that the polarization degree shows in general a random behavior, 
except in Cycle lib where it follows the 7?-band flux. Also, we found a lack of correlation between 
the polarization angle and the 7?-band flux. These results could be explained in terms of the 
presence of one or more variable polarization components. 

A possible scenar io that could explain our observational results is the "spine-sheath" model 



(Ghiselliniet al. 



2005[) . This model was proposed for TeV BL Lacs and invokes a slower external 
flow surrounding a fast spine. The slow layer (sheath) could be the result of the interaction of 
the jet walls with the ambient m edium. Another model that can be considered fo r our study 



is the two component model (see 



Havashida et al 



2008 



Krawczvnski et al 



2002). This model 



proposes the presence of two emission regions: a steady, quasi-stationary component, and another 
rapidly variable SSC component. Also, we want to mention that a multi-component scenario for 
lES 1959-1-650 has been propos ed to explain the observat i onal diff"erences found in the radio . 



optica l. X-ray, and y-ray bands (|Krawczynski et al 



2008h . 



2004 



Finer et al. 



2008 



Tagliaferri et al. 



4.1. Parameters of the stable component 

In terms of the normalized Stokes parameters, the absolute Stokes vector is defined as 
(Q, U)={q, u) I. To identify the p resence of a stable polarized component, we have used the method 



suggested by 



Jones et al. 



19851). In this work, the authors proposed that if the observed average 



values ((2), (f/)) in the absolute Stokes parameters plane Q-U deviate significantly from the 
origin, then a stable or a long-term (~ months-years) variable polarization component is present. 

For the case of a two-component model, we define the average values of Q and U as the 
stationary polarization component. These averages were calculated iteratively. First, we calculated 
the average of Q (or U) using all data, then the obtained value was recalculated again after 
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disca rding the outliers (> Scr). The iteration continued until no outliers remained (|lkejiri et al. 



20111). Figure |5] shows the Q-U plane for the stable component.The obtained average values of 
the absolute Stokes parameters (Q) and (U) are -0.10+0.01 mJy and 0.16+0.02 mJy, respectively. 
These average values appear as a black- filled star in the Q-U plane plot, and correspond to a 
stable component with constant polarization degree Pc = 4-.1 + 0.5% and polarization angle 
©c =151°+13°. The constant polarization degree has a dispersion crp^ =2.1%. 

4.2. Parameters of the variable component 

To find the origin of the variability behavior in Cycle II, we used the method proposed by 



Hagen-Thom et al. 



(|2008l) . Therefore, it is assumed that the variability within certain time interval 
is due to a single variable component. If the variability is caused only by its flux variations, while 
the relative Stokes parameter q and u remain unchanged, then in the space of the absolute Stokes 
parameters {/, Q, U} the observational points must lie on straight lines. The slopes of these lines 
are the relative Stokes parameters of the variable component. 

To estimate the variable component parameters, we looked for a possible linear relation 
between Q\s I and t/ vs / for the three relevant Cycles. For Cycles I and III no linear correlation 
between these parameters was found; rather, they appear to be randomly related. In contrast, 
for Cycles Ila and lib (see Figure |6l), our data show a linear tendency between these parameters. 
We made a least square fit to the data in order to find the slopes and the linear correlation 
coefficients r^ and r^,. The correlation coefficients obtained for these parameters in Cycle Ila are 
r^=0.913 and r„= -0.902, and the slopes m,^=0.228+0.028, m„=-0.408+0.054, respectively. The 
correlation coefficients for Cycle lib are r^=0.872 and r„=0.980, and slopes mg=0. 304+0.044, 
m„=0. 275+0.015, respectively. The Stokes parameters for the variable components qvar and 
Uvar found are given in Table |3] In this table, columns (2) to (7) show the parameters q^ar, rq, 
Uvar, fu, Pvar ^nd 6var, respectively. It is important to note, that Q vs I appear correlated during 
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Cycle Ila and Cycle lib. On the other hand, U \s I are anti-correlated during Cycle Ila and 
correlated in Cycle lib. The polarization degree found for the variable component in Cycle Ila is 
Pvar = (46.8 ± 4.9)%, with a polarization angle 9yar = 150° + 4°. For Cycle lib, the polarization 
degree is p^ar = (41.0 + 3.4)% and the polarization angle is O^ar = 21° + 4°. 

5. DISCUSSION 

5.1. Two-component Model 

The Q-U plane built from our data shows that the average values are shifted from the origin 
(Figure lU. Therefore, we can infer the presence of a constant or stable component that we 
assume associated with the relativistic jet, and also a variable component that can be related to 
the propagation of the shock. The observed polarization would be the result of the overlap of 
these two opt ically thin synchrotron components, according to the following expressions (see 



Holmes et al. 



1984 ): 



2 Pams + PLK/c + '^PvarPconshlc COS 2^ 
(1 +Ivlc) 

and 

tan20= , (17) 

Peons COS /.tfcons "I" Pvar^v/c COS Z.Ui,ar 

where p is the observed polarization degree and 6 the observed position angle. Here, ^ = 9cons-dvar 
and /v/c is the flux ratio between the variable to the constant component. 

To solve these equations, we used the previously found parameters P^ and 0^ as the 
parameters for the constant component p^ons and Ocons, respectively. To determine the ratio for the 
variable to constant component flux (7^/^ = Ivar/hons), we used data from the night JD 2455685 
where the contribution to the polarization of the variable component is almost zero, and the 
polarization degree is 1.5+0.8%. In this way, we fitted the observed polarization finding that 
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the ratio /v/c= 1.79+0.24. This corresponds to a constant component flux value of /, = 1.47 ±0.19 
mJy. With this result we solved the equations directly and found the values of p^,ar, O^ar and 
Ivar for each night. These results are presented for Cycles I and II in Table H and are plotted in 
Figure |7l where the polarimetric parameters of the variable component are shown during the low 
and high activity phases. In this figure it can be seen that during the low-activity phase (Cycle 
I) the polarimetric parameters of the constant polarized component, that are represented with a 
horizontal dashed-line, dominate over the variable polarized one (black-filled points). Therefore, 
for this cycle the observed polarization degree represented by empty squares, is much lower 
than the polarization degree of the variable component with a maximum diff'erence of 42%. 
Meanwhile, for the polarization angle the corresponding values for the variable and observed 
component, the difference is ~ 8%. For Cycle II, we can see that the observed polarimetric values 
(polarization degree and angle) are very similar to the variable component. This result shows 
that the variable polarization component dominates over the constant component. Using data 
from Table 121 columns (7) to (9), we find that the polarization degree is greater and more variable 
during the low-activity phase. 

From this analysis we can infer that the observed polarimetric behavior in the i?-band can be 
interpreted as the superposition of two optically-thin synchrotron components, one stable and the 
other variable. The p olarization degree p ^^„, ~ 4% derived for the stable component is identical 

(|2004|) . while the optical position angle Ocons shows value that 



to the value found by 



Bondi et al 



coincides with the ob served 9 ~ 1 50°. This va lue is similar to the projected position angle of 



the radio jet found by 



Finer et al. 



( 20081 



2010() . On the other hand, the source shows a larger 



polarization degree during the low-activity phase in comparison with the high-activity phase. The 
maximum polarization angle vari ations during the Iq w-ac tivity phase are ~50 °. These variations 



are similar to the value found by 



Rector et al. 



(2003) and 



Bondi et al. 



(120041) for the projected 



opening angle of the radio jet, while the maximum variations found during the high-activity phase 
are 10° and 18°, for Cycle Ila and Cycle lib, respectively. 
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The agreement between our optical results and those at radio bands (see above) suggests 
that there is a common magnetic field component for both the optical and radio bands, that 
is persistent in time and whose direction is transverse to the flow, as expected from a tangled 
magnetic field compressed by a shock. The polarimetric c haracteristics found in t his work are 



within the fra mework of a spine - sheath model proposed by 



Ghisellini et al, 



this source by .Finer et al. 



(l2005h and used for 



(|2010|) . In this context, the stable component would be associated with 



the sheath and the variable component with the spine. 



5.2. Standing Shock Scenario 



A standing shock scenario has already been proposed by iTagliaferri et al.l (|2008() 



Havashida et al. 



(120081) to explain the variable behavior of this object. A stationary shock wave 
system in a relativistic jet can be caused by the pressure imbalance between the jet and the 
surrounding ambient medium. Oscillations of the jet width occur as the overpressured medium in 
the jet expands until its pressure falls below the ambient pressure. This underpressured plasma 
then contracts under the influence of the external medium until high pressure is restored. This 
leads to the formation of a system of standing o blique shock waves, perhaps terminating in 



a strong shock perpendicular to the jet flow (see 



Gomez et al 



1997 



Sokolov et al 



2004 and 



references therein). To identify the origin of the variations in flux and polarization, we must 
consider that the two flares observed in 2009 and 2010 had time-scales of ~months. For this 
object, we have assumed that the radiation originates in the same region of the jet in a standing 
shock. 

Let us assume that the variable i?-band flux component is produced by a shock. In the 
one-zone homogeneous SSC model, the source is assumed to be a spherical blob of radius r^ 
moving in a turbulent plasma, with a constant jet's Lorentz factor F^, such that the emission region 
size is considered stable. Thus, in the observer's reference frame, the flux of the shocked region is 



21 



amplified as: 

F = Fqv-"5^^*''^6'^'^^"^ , (18) 

where 6 = [ry(l - p cos <!>]"' is the jet's Doppler factor, /? = (1 - Y^^Y^^ its global velocity in units 
of speed of light, and O is the viewing angle. The factor 6' is the Doppler factor of the shocked 
plasma in the shock's front reference frame. To determine 6' it would be necessary to know the 
plasma's state equation, nevertheless it cannot be determined since most of its parameters are 
unknown. However, without loss of generality it can be assumed that the plasma speed in the 



shock's reference frame s 



Hagen-Thom et al. 



lould be much less than the speed of light (« c) , assuming that 6' ~ \ 



2008 ) 



The observed degree of polarization depends on the shock's viewing angle ^, the spectral 
index or, and the ratio of densities of the shocked region to the unshocked region rj = rishockhunshock 



(IHughes & Milled 



1991 



and 



»P = tan ' 



Qr + 5/3 2-(l-77-2)sin2»P 
sinO 



>oso- ^i-rf)^ 



(19) 



(20) 



Following 



Tagliaferri et al 



(|2008|) . we assumed a bulk Lorentz factor F,- =18 for th e SED of 



lES 1959+650. We also have used the value of Uox = 1-64, given by 



AbdoetaL 



(l2010ch 



From equation (|T8l ) we can estimate the Doppler factor as a function of time. The value 
of Fq is determined by Fq = F,„a.xy" /So*"\ where F^a,: is the maximum observed flux and do is 
obtained from Oq, which is calculated from equations (fT9l) and (|20|) for the maximurn value of the 



polarization degree of the variable component (see Table |3) . From 



Hughes & Millej (Il99lh we 



adopt T] = 2.3, and *P = n/l which is the minimum possible compression that produces a degree 
of linear polarization as high as 47%. This yields *Po = 72.9°, Oq = 2.35°, and 5o=23.3 at the 
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maximum polarization (see Table |5]). 

Applying equation (fTSi) to the photometric data, we estimate the Doppler factor 6{t) and thus, 
the viewing angle of the jet 0(0 as functions of time. The latter allows us to estimate the viewing 
angle of the shock *P(0 (eq. [201) . Figure [8] shows the values of the Doppler factor, viewing angle of 
the jet and the viewing angle of the shock derived from the observed flux variability in Cycle 11. 
Applying equation (fT9l ) to the polarimetric data, we obtain the value of the plasma compression 
as function of time. The value of 77(f) is shown in Figure [8] as well. All these parameters are 
presented in Table |6l where we give in column (1) the Julian Date and in columns (2) to (5) the 
values obtained for 6{t), 0(f), ^(t) and 77(f) along with their errors, respectively. 

From Figure [81 it can be seen that when the source shows its maximum brightness (14.08 
mag, JD 2455002), the Doppler factor reaches 23.3, while in the minimum (14.74 mag, JD 
2455150) it is 20.4. This corresponds to a maximum variation of A6 = 2.9. The viewing angle 
of the jet O is anti-correlated with the flux, showing a minimum value of 2.35° when there is 
a maximum in brightness, and a maximum value of 2.78°. Therefore, it shows a maximum 
variation of AO= 0.43°. These small variations can produce large flux variations while Yj remains 
constant. In the state of maximum brightness, the viewing angle of the shock *P ~ 73° undergoes 
its maximum aberration due to relativistic efi'ects, deviating A¥ ~ 17° with respect to the shock 
transverse plane, while in a state of minimum brightness the shock tends to align perpendicularly 
to the axis of the radio jet. This finding may explain the high degree of polarization found in the 
low activity phase. That is, the amplification of the magnetic field components parallel to the 
shock due to Doppler eff"ect is not enough to considerably increase the polarization degree. In this 
context, the shock could be afl^ecting an emission zone with a magnetic field almost parallel to the 
jet axis. Another possibility is that the polarization degree may depend more on other factors, for 
example, the degree of compression in the electron density or depolarization by other components. 
We find the maximum compression of the plasma 77=1.145 when the polarization degree reaches 
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its maximum value of 10.3% (JD 2455062). We got the minimum compression factor ?7„,;„= 1.052 
when the polarization degree has a minimum value of 3.8% (JD 2454974). Thus, small changes 
in the compression factor At; « 0.093 can produce large changes in the polarization degree. 

5.3. Magnetic Field Structure 



The lifetime of the synchrotron electrons for a given frequency v in GHz (|Hagen-Thom et al 



20081) is: 



.2/ 1+^ '"" 



r,,, = 4.75 X 10^ d . (21) 

From our data we get t^.ar ~ t,nin = 9.74 ± 1.17 d, and Sq = 23.3. From here, we obtain the 
magnetic field intensity B = 0.061+0.005 Gauss, and an upper limit to the emission region size of 
rh < ct,ni,A)l{l +z) = (5.61+0.68)xlO" cm. 

The behavior of the variable polarization strongly depends on the amount of ordering of 
the magnetic field inside the source. The polarimetric properties observed in inhomogeneous 
sources at a given time and frequency will result from the integrated characteristics of all the 
different emitting regions. These inhomogeneities produce a depolarization degree in the emission 
source, revealing possible asymm etries in the magnetic field structure. In this work, we used the 



model proposed by 



Jones 



(|1988|) . where the inhomogeneities are produced by hydromagnetic 
turbulence, breaking the coherence beyond a characteristic distance Ig. Inside the emitting region, 
the magnetic field is modeled as a net of cubic cells with size Zg, which are subsequently divided 
into cells with size /b/2, until the observed results are fitted. Each ^-cell contains magnetic field 
vector's Bk, so the field's structure is considered to be frozen inside the plasma as it moves through 
the jet. Except for the influence of shocks, the field turbulence is isotropic in the fluid frame. The 
total magnetic field is the sum of all the field vectors inside the net. In the shock, both the density 
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and the field lines are compressed in such a way that the magnetic field components parallel to 
the shock's plane are amplified by r elativistic beaming, whereas the normal components are less 



influenced by relativistic eff"ects (see 



Hughes et al 



1985 ). 



Most of our data show that there is a lack of correlation between the total flux and the 
polarization degree. This result suggests that the evolution of the magnetic field is decoupled 
from the acceleration of the particles produced by the shock. An alternative way to obtain the 
structure of the magnetic field is through the calculation of the coherence length of the magnetic 
field at large scale Ib = (KUo/cTpy^^^l, where Ho = 0.75 is the fraction of polarization in a region 
with magnetic field perfectly ordered, k the degree of intrinsic pol arization of the source, cr^ the 



Jones et al. 



19851: 



Jones 



constant polarization dispersion, and I the emission region size (see 

Assuming that k is the maximum value of polarization for the variable component (/c=47%), 

taking cTp = 2% and / ~ r^, we obtain Z^ = 1.31 + 0.16 x 10'^ cm. 



1988). 



Observational and theoretical studies in radio (|Hughes et al 



1985 



Marscher & Gear 



1985h . 



show that variations in polarization on time-scales of several months are directly related to the 
distance along the jet traveled by the relativistic shock in the time between two extremes of 
the polarization l ight curve. If r^ is this distance, then the observed time scale At is given by 



(|Oian et al. 



19911) 



At 



r.(l+z) 



where c/3s is the speed of the shock, 6^ and Y^ are the Doppler fact or and Lorentz fac t or of 



(22) 



the shock, r espec t ively. Using values derived in radio at 43 GHz by 



Finer & Edwards 



Finer et al. 



( 2008L 



(2004); 



2010|) for the shock speed (J3s = 0.1), Lorentz factor of the shock (F^ « 3), and 
Doppler factor of the shock (5, ^ 2Ys), and taking the minimum observed time scale between two 
extremes of the polarization light curve, Ar=29d (Cycle III), we obtain r^ « (1.29 + 0.26)xl0^^ cm 
for the distance traveled by the shock. This distance is consistent with the coherence length Ib, 
suggesting a connection between the monthly variations observed in the polarization degree and 
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the spatial changes in the magnetic field induced by the inhomogeneities in the jet. All these 
derived parameters are summarized in Table |5] 



From our po 
at radio bands by 



arimetric analysis we find that our results are consistent with the results found 



Finer et al. 



(|2010l) . This work suggests that the behavior of the polarimetric 
parameters in radio are in good agreement with the spine- sheath model when applied to 
1 ESI 959+650. In this context, the electric vector position angle (EVFA) shows two difl'erent 
positions, one parallel along the jet axis, and another orthogonal at the edges of the jet. These 
authors also find that the polarization degree decreases as the jet is aligned toward the observer 
and increases at the edges of the jet. 



6. SUMMARY 

During the optical polarimetric observations of lES 1959+650 the object showed in the 
i?-band light curves mainly a period of low activity phase (2008), and later two active phases 
(flares) in 2009 and 2010. The time elapsed between the two maxima (14.08 mag at JD 2455002 
and 14.28 mag at JD 2455505) was ~500 d. From a detailed photometric and polarimetric analysis 
we have found the following results: 

1 . The minimum variability timescale of the i?-band flux was found to be ~10 d for the 



entire peri od of observations. Thi s result is similar to the value found by 



Krawczvnski et al. 



(120041) and 



Tagliaferri et al. (2008|) . We also found a maximum brightness variation of AmR=1.12 
mag in Ar=329 d. During the first maximum (7?= 14.08 mag) the brightness increases at a rate 
of 0.026 mJy/d and decreases at a rate of 0.025 mJy/d. During the second maximum (i?= 14.28 
mag) the source displayed a steady brightness increase from minimum brightness of i?=15.2 
mag to maximum brightness of 7?= 14.28 mag, with a rate of 0.020 mJy/d. In the entire period 
of observations, the source showed a minimum and maximum brightness of i?=15.2 mag and 
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i?= 14.08 mag, respectively. 

2. In general, there could not be found any correlation between the polarization degree 
and the optical flux, except when the source started to decrease its brightness after the first 
maximum occurred. This can be seen in Figure [I] during Cycle lib (2009 August 18, JD 2455062 
to November 16, JD 2455152). Therefore, only for Cycle lib a correlation between the optical 
R-band flux and the degree of linear polarization was found, with coefficient of correlation 
rpoi = 0.984 + 0.025. The maximum polarization degree of 12.2% was observed during the low 
activity phase, and the maximum variation of 10.7% occurred on a timescale of 2.8 yr. It is worth 
to mention that the polarization degree shows a more variable behavior (// ~ 32%) with a larger 
amplitude (Y ~ 128%) during the low-activity phase than during the higher activity phases where 
ju ~ 19% and Y ~ 64%. The minimum polarization degree time-scale is ~29 d. Furthermore, 
the source showed a larger average polarization degree during the low-activity phase (6.9%) in 
comparison with the high-activity phase (5.2%). No polarization angle rotation is observed during 
the monitored period of time. Instead, the source presented a preferential average position angle 
of 153+16°. This preferred position angle of the optical polarization is in g ood agreement with the 



projected position angle of the parsec scale jet found by 



Finer et al. 



(1201 oh from 43 GHz images. 



The lager dispersion of ~50° for the polarization angle was found during the low-actiyity ph ase. 



Rector et al. 



(2003). In 



This value is similar to the projected opening angle of the jet reported by 
the high-activity phase the dispersion turned out to be lower, i.e. ~18°. 

3. From the analysis done with the Stokes parameters {Q, U, /}, there follows the existence of 
two components that contribute to the polarized flux, one stable and another variable. The stable 
component has a constant polarization degree pc„ns = 4% and a polarization angle 9rn„^ ~ 150° . 



Bondi et al. 



(2004). 



These results are consistent with radio interferometric observations done by 

In this context, the polarimetric behavior is consistent with a spine-sheath structure of the jet 

where the variable component can be associated with the central jet structure taking the form of 
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spine and the constant component with a stable structure in the form of sheath surrounding the jet. 
4. The 7? -band observed variations can be explai ned in the framewo rk of a standing- shock 



model. Using the spectral index U nr = 1.64 given b y 



factor of the jet r^ = 18 found by 



Tagliaferri et al. 



Abdo et al. 



( 2010d) . and a bulk Lorentz 



(l2008h . a Doppler factor 6o = 23.3, and a 



viewing angle of the jet Oq = 2.35° were estimated. From our minimum variability timescale 
tmin we obtained the magnetic field intensity associated to the jet B = 0.06 Gauss, and the size of 
the emission region r^ = 5.6 x 10'^ c m. The magnetic field int ensity obtained in this work is in 



agreement with the value reported by .Krawczynski et al. 



(|2004r) . although the size of the emission 



region estimated by these authors (r^ « 1 .4 x 10 cm) is lower by an order of magnitude. 

5. Temporal variations of some parameters of the relativistic jet of lES 1959-1-650 were 
derived: (a) using photometric data, the Doppler factor 6(t), the viewing angle of the shock *F(0, 
and the viewing angle of the jet 0(0; (b) using polarimetric data, the ratio between the shocked to 
unshocked densities rj{t) and the coherence length of the magnetic field /e= 1.3 1x10" cm. 



6. From our data, and using values derived by 



Finer & Edwards 



(l2004h : 



Finer et al. 



(12008 . 



20101) for the physical parameters of the shock yS^ and F.,, we obtained the distance traveled by the 
shock r^= 1.29x10" cm. This distance is consistent with the field turbulence scale Ib, suggesting 
a connection between variations of the polarization degree and the spatial changes in the magnetic 
field. 

In a future work, an optical color polarization variability study will be carried out with the 
aim of studying the behavior of polarization with frequency. This will enable us to establish 
whether changes in the polarization degree found in this work are due to inhomogeneities in the 
density of particles in the emission zone or to changes in the magnetic field direction in the jet. 
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Fig. 1. — Light curves showing that the object displayed variations with timescales of -days- 
months. In top panel, the 7?-band (mag); middle panel, the polarization degree P(%); and bottom 
panel, the polarization position angle 6{°) vs Julian Date. Data are from 2007 October 18 to 2011 
May 05. Each point has an associated error bar and corrected by the contribution of the host-galaxy, 
see text. 
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Fig. 2. — A blow-up of the polarimetric parameter shown in Figure [T]for Cycles Ila and lib. There 
can be appreciated a clear correlation between the i?-band flux and the polarization degree during 
Cycle lib, while the polarization angle varied around the average value of ~ 150°. 
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Fig. 3. — Correlation between the i?-band flux and polarization degree for Cycle lib (2009 August 
14 to November 18). The Pearson's correlation coefficient is rpo/=0.984±0.025. There is a positive 
correlation between both parameters. 
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Fig. 4. — Top Panel: Correlation between the polarization angle and the 7? -band flux. Bottom 
panel: Correlation between the polarization angle and the polarization degree. In both panels there 
can be noticed a preferential tendency of the position angle ~ 150°. 
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Fig. 5. — Q-U absolute Stokes parameters plane obtained using all data. The solid black- filled 
star marks the mean constant value and the existence of a stable polarization component. The 
obtained average values are: (Q)= - 0.10+0.01 mJy and (U) = 0.16+0.02 mJy. The polarization 
degree and position angle values obtained for the constant component are Pc = (4.1 + 0.5)% and 
0c = 151° + 13°, respectively. The constant polarization degree dispersion is crp=2.\%. 
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Fig. 6. — Left: linear correlation between the Stokes parameters Q vs I (top panel), and U vs 
I (bottom panel) for Cycle Ila (2009 Apr-Jun). The correlation coefficients obtained for these 
parameters are r^=0.913 and r„= -0.902, and the slopes m^=0.228±0.028, m„=-0.408+0.054, re- 
spectively. Right: linear correlation between Q vs I (top panel), and U vs I (bottom panel) for 
Cycle lib. The correlation coefficients for this Cycle lib are r^=0.872 and r„=0.980, and slopes 
m^=0.304+0.044, m„=0.275+0.015, respectively. 
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Fig. 7. — Data collected during the low activity state in 2008 (left) and later during the high activity 
state in 2009 (right). They were fitted using a two-component synchrotron model. Squares show 
the superposition of a constant (dashed lines) and a variable polarized component (solid dots) using 
the addition law for polarized sources. In both figures, top panel shows the polarization degree and 
bottom panel shows the polarization angle. The model values are presented in Table 4. The values 
for the constant component presented in the plots are pams = (4.1 ± 0.5)% and Bcons = 151° + 13°. 
During the high activity state in 2009 (right), the variable component dominates over the constant 
component. 
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Fig. 8. — Temporal variability of some physical parameters related to the relativistic jet kinematics 
of lES 1959+650 derived in this work. From top to bottom panels: the Doppler factor, 6(ty, the 
angle between the jet axis and the line of sight, 0(0; viewing angle of the shock, ^(t); and the 
compression factor of the shocked to the unshocked plasma, ri(t), during the high activity state in 
2009. 



43 



Table 1 . POLARIZATION AND PHOTOMETRY IN THE 7?-BAND FOR lES 1959+650 

JD p ep e ee R 6r flux e/iux 

2450000.00+ (%) (%) (°) (°) (mag) (mag) (mJy) (mJy) 

4392.6685 8.5 0.6 147 02 14.71 0.02 4.04 0.10 

4396.6392 4.1 0.6 179 04 14.73 0.02 3.95 0.06 

4397.6553 5.6 0.7 167 03 14.66 0.02 4.22 0.06 

Cycle I 4591.9907 5.9 0.9 146 03 14.60 

4593.9912 6.0 1.0 138 03 14.77 

4594.9883 6.6 0.8 135 03 14.71 

4620.9575 3.9 0.8 144 04 14.69 

4621.9854 4.5 1.3 143 06 14.79 

4622.9287 3.3 0.9 160 07 14.71 

4623.9331 4.4 0.9 159 06 14.69 

4624.9839 7.4 1.0 164 03 14.79 

4655.9331 9.7 0.9 162 02 14.74 

4656.9160 12.2 0.9 142 02 14.78 

4660.9390 6.9 0.7 151 02 14.73 

4661.8857 6.5 0.8 153 03 14.77 

4681.9053 8.3 0.8 181 03 14.72 

4683.8711 4.4 0.8 151 04 14.73 

4707.7686 9.1 1.0 132 04 14.71 

4712.7568 6.2 0.9 137 05 14.60 

4737.6260 4.9 0.7 154 03 14.68 



0.02 


4.46 


0.07 


0.03 


3.82 


0.10 


0.02 


4.04 


0.09 


0.02 


4.09 


0.09 


0.03 


3.74 


0.10 


0.01 


4.04 


0.06 


0.01 


4.09 


0.06 


0.02 


3.74 


0.06 


0.02 


3.91 


0.06 


0.02 


3.78 


0.06 


0.01 


3.95 


0.06 


0.01 


3.82 


0.05 


0.01 


4.00 


0.06 


0.01 


3.95 


0.06 


0.02 


4.04 


0.10 


0.02 


4.46 


0.07 


0.01 


4.13 


0.06 



44 



Table 1 — Continued 



JD p ep e €e R e« flux efiux 

2450000.00+ (%) (%) (°) (°) (mag) (mag) (mJy) (mJy) 



4738.6494 5.9 0.7 168 03 14.66 

4764.6245 8.8 0.8 142 03 14.64 

4766.6401 6.6 0.8 148 04 14.63 

4767.6763 8.1 1.0 154 04 14.63 

4768.6421 6.4 0.8 142 04 14.68 

4771.6611 10.1 0.8 126 03 14.77 

4773.6025 10.3 0.9 186 03 14.75 

4774.6011 7.7 0.7 139 04 14.83 

4803.6006 5.1 1.0 153 05 14.83 

4804.5962 5.9 1.6 137 06 14.91 

Cycle Ila 4944.9717 7.3 0.5 144 03 14.34 

4973.9253 5.8 0.8 143 04 14.31 

4974.9087 3.8 0.8 145 05 14.26 

4975.9033 4.1 0.7 138 04 14.28 

4976.8965 4.1 0.7 142 03 14.31 

4977.9121 4.9 0.3 145 02 14.33 

4979.8735 4.9 0.4 146 02 14.40 

5002.9106 6.5 0.4 157 02 14.08 

5003.8867 5.7 0.4 153 02 14.11 

5004.9136 5.8 0.4 153 02 14.11 



0.01 


4.22 


0.06 


0.02 


4.27 


0.09 


0.01 


4.32 


0.06 


0.02 


4.32 


0.07 


0.01 


4.13 


0.06 


0.01 


3.82 


0.05 


0.01 


3.86 


0.06 


0.01 


3.61 


0.05 


0.01 


3.61 


0.05 


0.02 


3.33 


0.06 


0.01 


5.65 


0.08 


0.01 


5.84 


0.08 


0.01 


6.09 


0.08 


0.01 


5.96 


0.08 


0.01 


5.84 


0.08 


0.01 


5.71 


0.07 


0.02 


5.36 


0.12 


0.02 


7.19 


0.17 


0.02 


6.97 


0.16 


0.02 


6.97 


0.16 
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Table 1 — Continued 



JD p €p e 60 R eR flux 6fiux 

2450000.00+ (%) (%) (°) (°) (mag) (mag) (mJy) (mJy) 

5005.8696 5.6 0.4 151 02 14.16 0.02 6.69 0.16 

5006.8984 5.7 0.5 159 02 14.24 0.03 6.22 0.17 

5007.8599 4.9 0.4 151 02 14.17 0.03 6.62 0.17 

5008.8726 4.3 0.4 154 03 14.11 0.02 6.97 0.16 

Cycle lib 5059.7646 10.1 0.5 146 02 14.44 

5060.7266 9.4 0.5 147 01 14.42 

5062.7378 10.3 0.5 154 01 14.44 

5063.7241 9.4 0.5 147 02 14.43 

5064.7183 9.3 0.4 148 01 14.46 

5092.6494 7.6 0.6 139 02 14.49 

5093.6284 7.7 0.6 139 02 14.48 

5094.6235 7.4 0.6 138 02 14.56 

5095.6250 6.7 0.6 136 02 14.55 

5096.6221 6.8 0.6 139 02 14.56 

5097.6221 7.2 0.6 145 02 14.54 

5098.6680 7.7 0.6 144 02 14.56 

5122.6343 6.4 0.6 155 02 14.60 

5123.6606 6.9 0.5 153 02 14.56 

5124.6235 7.4 0.6 164 02 14.50 

5150.6030 7.0 0.7 171 02 14.74 



0.02 


5.14 


0.12 


0.02 


5.25 


0.12 


0.02 


5.14 


0.12 


0.02 


5.19 


0.12 


0.02 


5.08 


0.11 


0.01 


4.92 


0.07 


0.01 


4.98 


0.07 


0.01 


4.61 


0.06 


0.01 


4.66 


0.06 


0.01 


4.61 


0.06 


0.01 


4.71 


0.06 


0.01 


4.61 


0.06 


0.02 


4.46 


0.10 


0.01 


4.61 


0.06 


0.01 


4.87 


0.07 


0.02 


3.91 


0.09 
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Table 1 — Continued 



JD p ep e ee R sr flux ejiux 

2450000.00+ (%) (%) (°) (°) (mag) (mag) (mJy) (mJy) 

5152.5942 5.3 0.6 167 02 14.72 0.01 4.00 0.05 

5153.5879 6.1 0.6 164 02 14.71 0.01 4.04 0.05 

5154.5879 5.9 0.6 164 02 14.72 0.01 4.00 0.05 

Cycle III 5325.9937 2.5 0.7 157 08 15.14 

5327.9604 4.1 0.7 146 05 15.13 

5328.9199 3.5 0.6 152 05 15.09 

5329.9141 3.8 0.7 149 05 15.09 

5330.9302 3.1 0.6 143 05 15.13 

5331.9224 3.0 0.6 148 05 15.20 

5360.8804 4.0 0.4 147 03 14.99 

5362.9268 5.8 0.3 116 02 14.99 

5363.9282 3.4 0.3 178 03 15.00 

5364.9272 8.9 0.4 151 02 15.04 

5387.8569 6.1 0.3 186 02 14.86 

5390.8325 3.8 0.3 154 02 14.74 

5391.8589 3.5 0.4 162 03 14.88 

5392.8271 3.2 0.4 150 04 14.84 

5393.8359 4.3 0.4 156 03 14.84 

5413.7930 5.2 0.4 164 02 14.75 

5414.7808 2.9 0.4 155 04 14.75 



0.02 


2.70 


0.06 


0.02 


2.73 


0.06 


0.02 


2.83 


0.07 


0.02 


2.83 


0.07 


0.02 


2.73 


0.06 


0.02 


2.57 


0.06 


0.02 


3.11 


0.07 


0.02 


3.11 


0.07 


0.02 


3.07 


0.06 


0.02 


2.97 


0.07 


0.02 


3.49 


0.08 


0.02 


3.91 


0.09 


0.02 


3.45 


0.08 


0.02 


3.57 


0.08 


0.01 


3.57 


0.05 


0.01 


3.86 


0.05 


0.01 


3.86 


0.05 
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Table 1 — Continued 



JD p ep e ee R Sr flux efiux 

2450000.00+ (%) (%) (°) (°) (mag) (mag) (mJy) (mJy) 



5415.7773 3.7 0.4 156 03 14.80 

5416.7695 4.7 0.4 161 02 14.78 

5417.7744 5.2 0.4 166 02 14.77 

5418.8096 4.6 0.4 171 02 14.80 

5419.8135 8.4 0.5 191 02 14.79 

5446.7095 5.4 0.5 151 02 14.72 

5447.6626 2.2 0.5 137 04 14.74 

5448.6694 1.5 0.8 155 10 14.69 

5449.6602 1.7 0.8 162 08 14.69 

5450.6621 2.2 0.7 153 07 14.72 

5452.6562 7.1 0.9 196 05 14.72 

5476.6138 5.2 0.6 182 04 14.64 

5477.6196 4.8 0.4 180 02 14.64 

5478.6240 6.6 0.6 162 03 14.63 

5479.6304 8.0 0.5 154 02 14.62 

5480.6123 2.2 0.4 203 04 14.61 

5504.5923 6.0 0.5 143 02 14.32 

5505.6133 3.2 0.5 160 03 14.28 

5506.6655 3.8 0.6 158 03 14.31 

5507.6108 4.4 0.5 162 02 14.32 



0.01 


3.69 


0.05 


0.01 


3.78 


0.05 


0.01 


3.82 


0.05 


0.02 


3.69 


0.08 


0.02 


3.74 


0.09 


0.01 


4.00 


0.05 


0.02 


3.91 


0.09 


0.01 


4.09 


0.06 


0.01 


4.09 


0.06 


0.02 


4.00 


0.09 


0.01 


4.00 


0.06 


0.02 


4.27 


0.09 


0.01 


4.27 


0.06 


0.01 


4.32 


0.06 


0.01 


4.37 


0.06 


0.01 


4.41 


0.06 


0.01 


5.78 


0.08 


0.01 


5.96 


0.08 


0.01 


5.84 


0.09 


0.01 


5.78 


0.08 
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Table 1 — Continued 



JD p ep e ee R sr flux 6fiux 

2450000.00+ (%) (%) (°) (°) (mag) (mag) (mJy) (mJy) 

5508.6191 3.6 0.5 170 03 14.29 

5656.0156 2.0 0.7 158 08 14.91 

5683.9736 2.8 0.6 124 04 14.74 0.02 3.91 

5684.9956 2.9 0.7 142 05 14.68 

5685.9941 1.5 0.8 101 12 14.69 

5687.9917 3.3 0.7 134 05 14.62 



0.01 


5.90 


0.08 


0.02 


3.33 


0.08 


0.02 


3.91 


0.09 


0.01 


4.13 


0.06 


0.02 


4.09 


0.09 


0.01 


4.37 


0.06 
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Table 2. VARIABILITY PARAMETERS FOR lES 1959+650 



Cycle Parameter Average Max Min A„ 



Y(%) ni%) r 



X^ 



(1) 



(2) 



(3) 



(4) 



(5) (6) 



(7) 



(8) (9) 



(10) 



All 


R(mag; 


) 14.64 + 0.24 


15.20 


14.08 


1.12 


- 


- 


- 


- 




F(mJy; 


) 4.38 + 1.02 


7.19 


2.57 


4.62 


105.25 


23.28 


0.47 


13652.44 




p(%; 


) 5.7 + 2.3 


12.2 


1.5 


10.7 


187.57 


39.87 


0.79 


1611.09 




e{°\ 


) 153 + 16 


203 


101 


102 


65.60 


10.26 


0.34 


3219.04 


I 


R(mag; 


) 14.72 + 0.07 


14.91 


14.60 


0.31 


- 


- 


- 


- 




F(mJy; 


) 3.98 + 0.27 


4.46 


3.33 


1.13 


28.14 


6.72 


0.14 


464.34 




p(%; 


) 6.9 + 2.2 


12.2 


3.3 


8.9 


127.58 


31.87 


0.58 


163.13 




0{°\ 


) 150+14 


186 


126 


60 


39.83 


9.42 


0.19 


526.78 


Ila 


R(mag 


) 14.23 + 0.10 


14.40 


14.08 


0.32 


- 


- 


- 


- 




F(mJy; 


) 6.29 + 0.60 


7.19 


5.36 


1.82 


28.64 


9.46 


0.15 


263.40 




p(%; 


) 5.2+1.0 


7.3 


3.8 


3.5 


61.49 


19.22 


0.32 


45.28 




e{°\ 


) 149 + 6 


159 


138 


21 


13.47 


4.17 


0.07 


73.07 


lib 


R(mag 


) 14.55 ±0.10 


14.74 


14.42 


0.32 


- 


- 


- 


- 




F(mJy; 


) 4.67 + 0.43 


5.25 


3.91 


1.34 


28.35 


9.24 


0.15 


598.68 




p(%; 


) 7.6+1.4 


10.3 


5.3 


5.0 


63.52 


18.78 


0.32 


138.64 




e{°\ 


) 151 + 11 


171 


136 


35 


23.10 


7.29 


0.11 


563.81 


III 


R(mag 


) 14.77 + 0.23 


15.20 


14.28 


0.92 


- 


- 


- 


- 




F(mJy; 


) 3.90 + 0.88 


5.96 


2.57 


3.39 


86.83 


22.52 


0.40 


6005.49 




p(%; 


) 4.1 + 1.8 


8.9 


1.5 


7.4 


176.22 


43.51 


0.72 


553.05 




e{°\ 


) 157 + 19 


203 


101 


102 


64.01 


12.39 


0.34 


1657.01 
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Table 2 — Continued 



Cycle Parameter Average Max Min A^^x Y{%) fj.{%) T x^ 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



Note. — There are no statistics 7, //, T and x^ for the magnitude due to its 
logarithmic character. 



Table 3 . POLARIZATION FOR THE VARIABLE COMPONENT OF lES 1959+650 IN 

CYCLES Ila and lib. 

Cycle q^ar r^ U^ar ru Pvar(%) OvarH 

(1) (2) (3) (4) (5) (6) (7) 

Ila 0.228 + 0.028 0.913 -0.408 + 0.054 -0.902 46.8 + 4.9 150 + 4 
lib 0.304 + 0.044 0.872 0.275 ±0.015 0.980 41.0 + 3.4 21+4 
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Table 4. THE VARIABLE COMPONENT FOR TWO-COMPONENT MODEL 



JD 


Pvar 


"var 


'^var 


2450000.00+ 


(%) 


(°) 


(mJy) 


(1) 


(2) 


(3) 


(4) 


4591.9907 


6.8 + 0.7 ] 


[45 + 05 


3.00 + 0.26 


4593.9912 


7.6+1.1 ] 


134 + 06 


2.35 + 0.29 


4594.9883 


8.4+1.1 ] 


[31+05 


2.57 + 0.29 


4620.9575 


3.8 + 0.5 ] 


[40 + 10 


2.62 + 0.29 


4621.9854 


4.9 + 0.7 ] 


139 + 11 


2.27 + 0.29 


4622.9287 


3.0 + 0.3 ] 


[67 + 14 


2.57 + 0.25 


4623.9331 


4.7 + 0.5 ] 


[63 + 10 


2.62 + 0.25 


4624.9839 


9.8 + 1.3 ] 


[67 + 05 


2.27 + 0.25 


4655.9331 


13.4+1.6 ] 


[64 + 03 


2.44 + 0.25 


4656.9160 


17.5 + 2.2 ] 


[41 + 03 


2.31+0.25 


4660.9390 


8.6 + 1.0 ] 


[51 +04 


2.48 + 0.25 


4661.8857 


8.0+1.0 ] 


[54 + 06 


2.35 + 0.25 


4681.9053 


12.2 + 1.4 ] 


[86 + 04 


2.53 ± 0.25 


4683.8711 


4.6 + 0.5 ] 


[51 +09 


2.48 + 0.25 


4707.7686 


12.5 + 1.7 ] 


[29 + 05 


2.57 ± 0.29 


4712.7568 


7.5 + 0.8 ] 


[33 + 07 


3.00 + 0.26 


4737.6260 


5.4 + 0.6 ] 


[55 + 07 


2.66 + 0.25 


4738.6494 


7.4 + 0.8 ] 


[73 + 05 


2.76 + 0.26 


4764.6245 


11.5 + 1.4 ] 


[40 + 04 


2.80 ± 0.29 
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Table 4 — Continued 



JD 


Pvar 


"var 


^var 


2450000.00+ 


(%) 


(°) 


(mJy) 


(1) 


(2) 


(3) 


(4) 


4766.6401 


7.9 + 0.8 


147 + 06 


2.85 + 0.26 


4767.6763 


10.2+ 1.1 


155 + 05 


2.85 + 0.26 


4768.6421 


7.7 + 0.9 


140 + 06 


2.66 + 0.25 


4771.6611 


15.0+1.9 


122 + 04 


2.35 + 0.25 


4773.6025 


15.9 + 2.0 


190 + 03 


2.40 + 0.25 


4774.6011 


10.5 + 1.4 


136 + 06 


2.14 + 0.24 


4803.6006 


5.8 + 0.8 


154 + 09 


2.14 + 0.25 


4804.5962 


7.8 + 1.3 


132 + 09 


1.87 + 0.26 


4944.9717 


8.5 + 0.6 


143 + 04 


4.19 ±0.27 


4973.9253 


6.4 + 0.5 


141 +05 


4.37 ± 0.28 


4974.9087 


3.7 + 0.3 


143 + 08 


4.62 + 0.28 


4975.9033 


4.2 + 0.3 


134 + 06 


4.49 + 0.28 


4976.8965 


4.2 + 0.3 


139 + 06 


4.37 + 0.27 


4977.9121 


5.2 + 0.4 


143 + 04 


4.25 + 0.27 


4979.8735 


5.2 + 0.5 


145 + 04 


3.89 + 0.31 


5002.9106 


7.2 + 0.5 


158 + 03 


5.72 + 0.36 


5003.8867 


6.1+0.4 


153 + 03 


5.50 + 0.36 


5004.9136 


6.2 + 0.5 


153 + 03 


5.50 + 0.36 


5005.8696 


6.0 + 0.4 


151+03 


5.22 + 0.35 



53 



Table 4 — Continued 



JD 


Pvar 


"var 


^var 


2450000.00+ 


(%) 


(°) 


(mJy) 


(1) 


(2) 


(3) 


(4) 


5006.8984 


6.3 + 0.5 


161+03 


4.75 + 0.36 


5007.8599 


5.1 +0.4 


151 +04 


5.15 + 0.36 


5008.8726 


4.3 + 0.3 


155 + 04 


5.50 + 0.36 


5059.7646 


12.5 + 1.2 


145 + 03 


3.67 + 0.31 


5060.7266 


11.5 + 1.1 


146 + 02 


3.78 + 0.31 


5062.7378 


12.7 + 1.2 


154 + 02 


3.67 + 0.31 


5063.7241 


11.4+ 1.1 


146 + 03 


3.73 + 0.31 


5064.7183 


11.4+ 1.1 


148 + 02 


3.62 + 0.31 


5092.6494 


9.3+0.8 


137 + 03 


3.45 + 0.26 


5093.6284 


9.3 + 0.8 


137 + 03 


3.51 ±0.26 


5094.6235 


9.2 + 0.9 


135 + 03 


3.14 + 0.26 


5095.6250 


8.2 + 0.8 


133 + 04 


3.19 + 0.26 


5096.6221 


8.3 + 0.8 


136 + 04 


3.14 + 0.26 


5097.6221 


8.6 + 0.8 


144 + 03 


3.25 + 0.26 


5098.6680 


9.5 + 0.9 


143 + 03 


3.14 + 0.26 


5122.6343 


7.6 + 0.8 


156 + 04 


3.00 + 0.29 


5123.6606 


8.3+0.8 


154 + 04 


3.14 + 0.26 


5124.6235 


9.0 + 0.8 


166 + 03 


3.40 + 0.26 


5150.6030 


9.6 + 1.3 


176 + 03 


2.44 + 0.28 
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Table 4 — Continued 



J^ Pvar "var ^var 

2450000.00+ (%) (°) (mJy) 

(1) (2) (3) (4) 

5152.5942 6.5 + 0.7 173 + 04 2.53 + 0.25 

5153.5879 7.5 ± 0.8 168 + 04 2.57 + 0.25 

5154.5879 7.2 ± 0.8 168 + 04 2.53+0.25 
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Table 5. PHYSICAL PARAMETERS FOR lES 1959+650 



Parameter 


Value 


Units 


B 


0.06 + 0.01 


Gauss 


So 


23.3 + 0.4 




^0 


72.9 + 6.1 


degree 


<Do 


2.35 + 0.26 


degree 


n 


5.61 +0.68 


10" cm 


rs 


1.29 + 0.26 


10" cm 


Ib 


1.31 +0.16 


10" cm 


i-min 


9.74 ± 1.17 


days 



Note. — The parameters were 
calculated using values derived by 



iner & Edwards 



( 2008 



(2004) and 



Piner et al. 



20101) for the shock speed Ps = 
0.1, and the Lorentz factor of the shock 
Ys ~3. In addition, we use the value for 
the optical spectral index a„y = 1 .64, re- 



ported by 



AbdoetaL 



(l2010cl) : and the 



value of the bulk Lorentz factor of the 
jetrj=18 from 



Tagliaferri et al 



(2008) 



[°] 


[°] 


(3) 


(4) 
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Table 6. VARIABLE PHYSICAL PARAMETERS FROM POLARIMETRIC DATA 

OBTAINED FOR lES 1959+650 

JD 6(t) 0(0 T(0 77(0 
2450000+ 

(1) (2) (3) (4) (5) 

4944.9717 22.1+0.4 2.52 + 0.06 76.8 + 1.3 1.102 + 0.008 

4973.9253 22.3 + 0.4 2.50 + 0.06 76.3 + 1.3 1.080 + 0.011 

4974.9087 22.5 + 0.4 2.47 + 0.06 75.6 + 1.3 1.052 + 0.012 

4975.9033 22.4 + 0.4 2.49 + 0.06 76.0+1.3 1.056 + 0.011 

4976.8965 22.3 + 0.4 2.50 + 0.06 76.3 + 1.3 1.056 + 0.010 

4977.9121 22.2 + 0.4 2.52 + 0.06 76.6 + 1.3 1.067 + 0.005 

4979.8735 21.9 + 0.4 2.56 + 0.06 77.6 + 1.3 1.067 + 0.006 

5002.9106 23.3 + 0.4 2.35 + 0.06 72.9 + 1.4 1.095 + 0.007 

5003.8867 23.1+0.4 2.38 + 0.06 73.4 + 1.4 1.081+0.006 

5004.9136 23.1+0.4 2.38 + 0.06 73.4 + 1.4 1.083 + 0.007 

5005.8696 22.9 + 0.4 2.41+0.06 74.1 + 1.4 1.079 + 0.006 

5006.8984 22.6 + 0.4 2.46 + 0.06 75.3 + 1.3 1.080 + 0.007 

5007.8599 22.9 + 0.4 2.41+0.06 74.3 + 1.4 1.068 + 0.006 

5008.8726 23.1+0.4 2.38 + 0.06 73.4+1.4 1.060 + 0.006 

5059.7646 21.7 + 0.4 2.59 + 0.06 78.3 + 1.3 1.142 + 0.008 

5060.7266 21.8 + 0.4 2.58 + 0.06 78.0+1.3 1.133 + 0.007 

5062.7378 21.7 + 0.4 2.59 + 0.06 78.3 + 1.3 1.145 + 0.007 

5063.7241 21.7 + 0.4 2.58 + 0.06 78.1 + 1.3 1.132 + 0.007 

5064.7183 21.6 + 0.4 2.60 + 0.06 78.4+1.2 1.130 + 0.007 



57 



Table 6 — Continued 



0(0 


¥(0 


[°] 


[°] 


(3) 


(4) 



JD 6(t) 0(0 W(t) r](t) 

2450000+ 

(1) (2) (3) (4) (5) 

5092.6494 21.5 + 0.4 2.62 + 0.06 78.9 + 1.2 1.105 + 0.008 

5093.6284 21.5 + 0.4 2.61+0.06 78.8 + 1.2 1.106 + 0.008 

5094.6235 21.2 + 0.4 2.67 + 0.06 79.9 + 1.2 1.101+0.008 

5095.6250 21.2 + 0.4 2.66 + 0.06 79.7 + 1.2 1.091+0.008 

5096.6221 21.2 + 0.4 2.67 + 0.06 79.9 + 1.2 1.093 + 0.008 

5097.6221 21.3 + 0.4 2.65 + 0.06 79.6 + 1.2 1.098 + 0.008 

5098.6680 21.2 + 0.4 2.67 + 0.06 79.9 + 1.2 1.105 + 0.008 

5122.6343 21.0 + 0.4 2.69 + 0.06 80.4+1.2 1.087 + 0.008 

5123.6606 21.2 + 0.4 2.67 + 0.06 79.9 + 1.2 1.094 + 0.007 

5124.6235 21.4 + 0.4 2.63 + 0.06 79.1 + 1.2 1.101+0.008 

5150.6030 20.4 + 0.4 2.78 + 0.06 82.3 + 1.2 1.095 + 0.009 

5152.5942 20.5 + 0.4 2.77 + 0.06 82.0+1.1 1.070 + 0.008 

5153.5879 20.6 + 0.4 2.76 + 0.06 81.8 + 1.1 1.081+0.008 

5154.5879 20.5 + 0.4 2.77 + 0.06 82.0+1.1 1.078 + 0.008 

Note. — Active cycle 2009 for lES 1959+650. Col. (1): Julian Day; 
Col. (2): Doppler factor; Col. (3): Angle between the line of sight and jet 
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axis; Col. (4): Viewing angle of the shock; Col. (5): Ratio of the density 
in the shocked region with respect to unshocked region. 



